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We show numerically that ultrashort self-defocusing temporal solitons colliding with a weak pulsed probe in the 
near-IR can convert the probe to the mid-IR. A near-perfect conversion efficiency is possible for a high effective 
soliton order. The near-fR self-defocusing soliton can form in a quadratic nonlinear crystal (beta-barium borate) 
in the normal dispersion regime due to cascaded (phase-mismatched) second-harmonic generation, and the 
mid-IR converted wave is formed in the anomalous dispersion regime between A = 2.2 — 2.4 pm as a resonant 
dispersive wave. This process relies on non-degenerate four-wave mixing mediated by an effective negative 
cross-phase modulation term caused by cascaded soliton-probe sum-frequency generation. 

OCIS codes: (320.7110) Ultrafast nonlinear optics; (190.5530) Pulse propagation and temporal soli¬ 

tons; (320.2250) Femtosecond phenomena 
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The optical soliton is remarkably robust as it can both 
retain its shape despite dispersive or dissipative effects 
and survive wave collisions. Yet it is quite susceptible to 
perturbations as it can shed phase-matched resonant ra¬ 
diation to a so-called soliton-induced optical Cherenkov 
wave (a.k.a. dispersive wave) when perturbed by higher- 
order dispersion [l-5j. On the other hand, the soliton 
can also act as a potential barrier when colliding with 
a linear (i.e. dispersive) wave, creating the analogy to 
an optical equivalent of an ’’event horizon” iS, 0. Such 
a collision can be well understood by generalizing the 
Cherenkov phase-matching condition - degenerate four- 
wave mixing (FWM) - to non-degenerate FWM where 
a soliton interacts with two linear dispersive waves [§- 
D. This interaction is mediated by cross-phase mod¬ 
ulation (XPM): the collision between the soliton (w s ) 
and a linear ’’probe” wave (w p ) can become resonantly 
phase-matched to a new frequency (the ’’resonant” wave 
w r ) according to the FWM phase-matching condition. 
This frequency-converts the probe to the resonant wave, 
which - when completely depleting the probe - gives rise 
to the peculiar appearance of the probe reflecting on the 
soliton: when the probe group velocity is higher than the 
soliton, the resonant wave group velocity will be lower 
than the soliton, and therefore travel away from the soli- 
ton after formation. This frequency-conversion process 
therefore in time domain leads to an apparent reflection 
of the probe on the soliton. 

Soliton-probe collisions have been studied in fibers, 
which have a positive Kerr nonlinearity. How¬ 
ever, through cascaded (strongly phase-mismatched) 


quadratic nonlinear interactions an effective negative 
Kerr-like nonlinearity may be generated (in bulk this 
corresponds to a self-defocusing effect) [l3|. As a con¬ 
sequence soliton formation requires normal dispersion 
[13L 14] , and Cherenkov phase-matching naturally occurs 
towards the red side of the soliton spectrum (15j|. This 
allows for efficient near- to mid-IR conversion jig, E3- 

We here study the collision of near-IR probe waves 
and self-defocusing solitons to generate long-wavelength 
resonant waves. The self-defocusing soliton may reflect 
the probe wave when the probe XPM term is nega¬ 
tive, and this is possible if the sum-frequency generation 
(SFG) between soliton and probe is detuned sufficiently 
away from its phase-matching point to induce a negative 
cascaded XPM term. We find that complete conversion 
from the probe to the resonant wave requires interaction 
with a higher-order soliton, and demonstrate a wide tun- 
ability of the wavelength of the mid-IR resonant wave by 
varying the probe center wavelength. 

The general FWM phase-matching condition is Q 

k\]n r ) = k S ol(u] r ) “t“ J (ujp ) k so \ (^'p)] (1) 

where ku n (u>) describes the dispersion relation of the lin¬ 
ear wave (in bulk media simply determined, e.g., by the 
Sellmeier equation). fc so i = fciin(w s ) + (w-w s )/u SiSO i + q S oi 
is the soliton dispersion relation; its non-dispersive na¬ 
ture is reflected in the fact it is simply a wave packet 
with a constant group velocity u ff>so i. Its accumulated 
nonlinear phase g so i will cancel out for the J = +1 case 
that we will focus on here. The parameter J switches be¬ 
tween the degenerate case (J = 0, Cherenkov radiation) 
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Fig. 1. (a) The dispersion relations for BBO in frequency 
domain, reported in the soliton group-velocity frame. Black 
line: left-hand side of Eq. 0, colored lines: right-hand side 
of Eq. ([TJ), taking a soliton at \ s = 1.1 pm and a probe at 
X p = 1.65 pm. ”A” and ”N” denote regimes with anomalous 
and normal dispersion, respectively, (b) The group velocities 
of a linear wave (black) and the soliton at A s = 1.1 pm (blue). 

and the non-degenerate case (J = ±1, where the pres¬ 
ence of the probe at frequency uj p invokes the FWM res¬ 
onant phase-matching condition). In a BBO quadratic 
nonlinear crystal (/3-barium borate, BaB 204 ) the res¬ 
onant waves are phase-matched in the mid-IR beyond 
A = 2.0 pm, as seen from the dispersion relations in Fig. 
|T] for the main case considered here, namely a 1.65 pm 
probe colliding with a 1.1 pm soliton. 

The BBO crystal is assumed cut for type-I second- 
harmonic generation (SHG), where two o-polarized pho¬ 
tons at the fundamental wave (FW) frequency lo\ gen¬ 
erate a second-harmonic (SH) e-polarized photon at the 
frequency W 2 = 2wi. The numerics use the nonlinear 
wave equations in frequency domain [l8| model; the val¬ 
ues of the x*- 2 ) and \tensor components were chosen 
from [l9| , and the Raman effect is neglected as it is usu¬ 
ally considered weak in BBO. Here we pump in the o- 
wave and through phase-mismatched SHG to the e-wave 
a nonlinear phase shift accumulates on the pump pulse, 
which we exploit to excite a self-defocusing soliton. Im¬ 
portantly, since we also pump with a weak o-polarized 
probe wave, the numerical model also includes any pos¬ 
sible interaction, such as sum- and difference fre¬ 
quency generation (SFG and DFG), and both inter- and 
intra-polarization (i.e. type 0, I and II) interactions. For 
simplicity the BBO mid-IR material loss is neglected. 

The collision is modelled by launching 
two co-propagating o-polarized fields -E 0 ,in = 
Ei t i n cos(wif) sech(t/Ti) + I3p i in cos (cjpt) sech[(t — t)/T p \, 
where r is the delay time between them and we will 
use identical input pulse durations T\ = T p . A 
self-defocusing soliton can form at uj s = oj\ if the 
following criteria are fulfilled (see [20] for more de¬ 
tails): (a) The effective Kerr self-phase modulation 
(SPM) nonlinearity ri 2 ™( w i) = n pSc( w i) + n 2 ,Kerr 
must be negative, and this is controlled by 
n|^£ c ( w i) oc — dl s /Ak^ G by making the SHG 
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Fig. 2. Evolution of the o-polarized field in (a) time domain 
(normalized intensity of the electrical field envelope) and (b) 
wavelength domain, (c) The dispersion curves of Fig. [lj a) 
in wavelength domain for A s = 1.1 pm (thick) and A s = 
1.075 pm (thin). Input pulses: 50 fs FWHM@Ai = 1.1 pm, 
J l in = 200 GW/cm 2 (soliton) and 50 fs FWHM@A P = 
1.65 pm, Ipp n = 5 GW/cm 2 , r = 150 fs (probe). The BBO 
crystal had A k^ G = 60 mm -1 (6 = 18.8° and ip = —90°). 

phase mismatch A k^ G = k e (2w\,6) — 2k 0 (wi) small 
enough; (b) the group-velocity dispersion (GVD) 
must be normal [k y 0 ; (wi) > 0], which in BBO 
means Ai < 1.488 pm; (c) the effective soliton order 
[SI N eS > 1, where iV e 2 ff = L D wi/ in |n!^(wi)|/c, 

h,in = £oni(wi)c|£;i i in| 2 /2 and L D = T%/k < £‘\ op). 

The soliton will reflect/scatter the probe if the soliton- 
induced XPM potential is a barrier. As the probe is 
launched in the anomalous dispersion regime, this re¬ 
quires a negative probe XPM effective nonlinear in¬ 
dex. This has the contributions j2ti}: n*™ = + 

!casc( w s + 9 ) + - Wp, 0), viz. the material 

Kerr XPM (which in BBO is identical to ri 2 ™ rr , as both 
waves are o-polarized), as well as type I cascaded soliton- 
probe SFG and DFG. These are similar in form to nf^ c 
[2p| and in our case we find that nJ/^ c (w. s — uj p ,6) is 
negligible. Thus, the material Kerr XPM effect and the 
cascaded SFG effect are the main contributions to the 
XPM sign and magnitude. Only in certain regimes may 
n % ; ™ < 0 [2C], Fig. 6], and the potential may even flip 
sign to become a hole. 

Fig.[2]shows the results from a typical simulation. The 
BBO crystal angle is suitably chosen to give nonreso¬ 
nant [231 negative SPM and XPM nonlinearities (which 
occurs between 17.5° < 6 < 20.0° 0, Fig. 6]). The 
soliton input intensity is chosen so N e g = 2.0, allowing 
a higher-order self-defocusing soliton to form. A weak 
probe is launched in the anomalous dispersion regime, 
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which from Fig. [lfb) implies that its group velocity is 
larger than the soliton. It is therefore suitably delayed at 
the input so the interaction occurs over realistic crystal 
lengths. The time plot in (a) shows the probe colliding 
with the trailing edge of the strong soliton at around 
10 mm. After the collision a reflected wave emerges; 
this is the resonant wave phase-matched to the soliton 
through the negative XPM nonlinearity. According to 
Fig. QJb) the resonant wave will have a lower group ve¬ 
locity than the soliton and this explains why it is trav¬ 
eling away from the soliton trailing edge. In wavelength 
domain (b) the normalized spectral density (SD, calcu¬ 
lated as S( A) = | A(u = 2ttc/X)\ 2 X 2 /c) shows that the 
probe is almost completely converted to the resonant 
wave, and this occurs between 10 and 20 mm propaga¬ 
tion. It is exactly in this propagation range the collision 
takes place in time domain. There is a good agreement 
between the predicted phase-matching frequency of the 
resonant wave, which is evident from the wavelength- 
domain phase-matching curves plotted in (c). 

We note that the soliton blue-shifts slightly during 
propagation as a consequence of cascading-induced self- 
steepening. This means that at the final stage uj s > wi, 
and this leads to a new set of soliton curves for the phase¬ 
matching conditions, indicated as thin lines in (c). As 
these represent solitons they are in frequency domain 
still straight curves, but they are now tilted instead of 
flat as the group velocity is different from that at u >±; in 
wavelength domain this means that they are no longer 
represented as straight curves, but this is simply due to 
the A oc 1/w relation. These curves explain how the res¬ 
onant wave is found slightly more red-shifted than the 
w s = u>i case predicted. We also see that the Cherenkov 
(J = 0) case can be seen in the spectrum, accurately pre¬ 
dicted by the blue-shifted soliton phase-matching condi¬ 
tion. This is the degenerate case, where the soliton alone 
becomes phase-matched to a resonant wave. 

Detailed spectral and temporal cuts are shown in Fig. 
EH (a) shows the SD at input, during (z = 15 mm) and 
after collision (z = 30 mm). Since the soliton order is 
above unity, the o-polarized spectra (thick lines) show 
that the soliton at collision is considerably extended to¬ 
wards the probe spectrum, and the final spectrum shows 
that the probe is almost completely depleted leaving 
only the J = +1 and J = 0 resonant waves. The o- 
polarized third harmonic is also evident, and in the e- 
polarized spectra (thin lines) the various SHG and SFG 
components are evident as well. The same cuts are shown 
in time domain focusing in (b) on the probe (using a 
band-pass filter) and in (c) on the resonant wave (us¬ 
ing a long-pass filter). The probe at 15 mm is around 
half depleted, giving most of its depleted energy to the 
resonant wave that is located at the same temporal po¬ 
sition. After 30 mm the weak probe does not show on 
a linear scale (in the plot it is amplified 10 times). The 
resonant wave is now delayed 400 fs and it is reduced in 
intensity and increased in time due to dispersion. It has 
a Gaussian profile since it is a linear and not a soliton 





Fig. 3. Cuts from Fig. [2] (a) the normalized o- (thick) and 
e-polarized (thin) SD; (b) the band-pass filtered probe and 
(c) the long-pass filtered resonant wave envelopes vs. time, 
(d+e) Normalized energy vs. z in the collision (full) and no 
collision (dashed) cases of (d) the soliton and (e) the entire 
e-polarized spectrum, the filtered probe and resonant waves. 

wave. Finally (d+e) show the energy, normalized to the 
total input energy, of the soliton, probe and resonant 
waves. The soliton initially looses around 2% of its en¬ 
ergy through SHG to the e-polarized SH [(e) also shows 
the total e-polarized energy], causing the initial ripples 
at z < 1 mm. The soliton-probe interaction occurs be¬ 
tween z = 5 — 20 mm, and the resonant wave builds up 
in energy. After 20 mm the probe is depleted. The en¬ 
ergy ratio (conversion efficiency) of the resonant wave to 
the probe is around 0.72, close to the limit posed by the 
photon-to-photon ratio ui r /ui p = 0.73 as dictated by the 
Manley-Rowe relation. The energies from a simulation 
where the probe never collides with the soliton (dashed 
lines) show as expected no energy at the resonant wave 
and the probe remains unaffected. 

Fig. U shows the probe and resonant wave energies vs. 
N e g (controlled by /yin) for (a) X p = 1.65 fim. fixed and 
four different / Pi i„ values, and (b) / Pi i„ = 10 GW/cm 2 
fixed and three different A p values. In (a) the complete 
depletion of the probe happens for intensities up to 
30 GW/cm 2 , all ending up at the same plateau, whose 
level is dictated by the Manley-Rowe relation (as oj p is 
fixed). For I p g n = 40 GW/cm 2 the probe conversion is 
incomplete, which is a result of nonlinear spectral broad¬ 
ening occurring before collision. It is therefore impor¬ 
tant that the probe does not experience any nonlinear 
phase shifts before the collision. In (b) a fixed moderate 
probe intensity was therefore used to ensure complete 
depletion of the probe. Since different probe frequen¬ 
cies are used the plateau levels vary, in accordance with 
the Manley-Rowe relation. Clearly probe depletion re¬ 
quires JV e g > 1 . Generally the resonant wave growth 
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Fig. 4. The probe and resonant wave energies (normalized 
to probe input energy) vs. at z = 40 mm. (a) Fixing 
A p = 1.65 /im and (b) fixing I P} i n = 10 GW/cm 2 . 


follows a logistic sigmoid function, whose slope scales as 
J~ in 5 and the midpoint N e g^ increases linearly with I p , j n 
and decreases linearly with X p . This latter scaling comes 
from the fact that as X p increases it approaches the zero 
group-velocity mismatch (GVM) wavelength, where the 
probe and soliton have identical group velocities, see Fig. 
mb). The detuning from GVM has traditionally been 
kept small to get a large reflection, with the dilemma 
that the resonant wavelength is almost identical to the 
probe [see again Fig'- [Ub)] and the collision will only take 
place through very long interaction distances EJ- Still 
a decent conversion should be possible even with large 
detunings from the zero GVM point (see e.g. [HI), and 
we intend to investigate this further in another paper. 

The emission wavelength of the nrid-IR resonant wave 
is tunable from A = 2.2 — 2.4 /im through changing the 
probe wavelength, see Fig. [5] Note that the nrid-IR wave 
remains quite significant even when probe-depletion is 
incomplete and when IR losses are taken into account. 
Besides, this process can be generalized to nrid-IR trans¬ 
parent crystals like LiNbOa mm, where the phase¬ 
matching point will also lie further into the mid-IR. 

Concluding, in a quadratic nonlinear crystal multi¬ 
ple cascaded nonlinear effects allow firstly the forma¬ 
tion of a self-defocusing near-IR soliton (through cas¬ 
caded SHG generating a self-defocusing Kerr-like SPM 
term), which when colliding with a weak near-IR probe 



Fig. 5. Long-wavelength part of the o -polarized SD for various 
probe wavelengths shown after z = 20 mm (except X p = 
1.75 fj,m where z = 30 mm) with (thick) and without (thin) 
IR losses of BBO. The dashed lines show the input states. 
The soliton was the same as in Fig. [5] and I p = 10 GW/cm 2 . 


induces a cascaded SFG nonlinearity generating a self- 
defocusing Kerr-like XPM term. This allows a resonant 
mid-IR wave to become phase matched. A complete 
probe-resonant wave conversion is possible when collid¬ 
ing with a higher-order soliton, a case which has not 
been considered before in Kerr systems, and for a given 
soliton wavelength the resonant wavelength is tunable by 
varying the probe wavelength. Obtaining phase match¬ 
ing further into the mid-IR is possible in other crystals, 
and interestingly the system allows to change the XPM 
term sign to study the barrier vs. hole potential effect. 

X.L., M.B. and B.Z. acknowledge the support from the 
Danish Council for Independent Research (11-106702). 

Appendix A: Theoretical background 

Let us first stress that the code we use models the o- 
and e-polarized electrical fields and include all possible 
X ^ and x ^ interactions. However, an analytical un¬ 
derstanding of the processes behind the interaction is 
more convenient in the slowly varying envelope approxi¬ 
mation. We will now consider the interaction of a strong 
soliton envelope E s at the frequency and a weaker 
(but not weak) probe envelope E p at the frequency u p . 
Let us focus on the two equations for the soliton and 
probe and disregard dispersion for the moment. For the 
soliton 


UA. + 


cn a (uj s ) 


V E U j E U j s — U j 


-*Afc DPG , 


3w s 


8 n Q (w s )c 


I T 7”* p 

d 3 ) 


[XspmI^s f + 2Xxpm\Ew p f ] — 0 


(Al) 


where n 0 is the o polarized refractive index, modelled by 
the BBO Sellmeier equations reported in [24|. Here we 
have only included the type I oo —> e SHG, SFG and 
DFG possibilities, all having the same effective nonlin¬ 
earity 


d e ff (9,<j)) =d 3 i sin 9 — d 2 2 cos 9 sin 3<j> (A2) 

For the cut we use this is optimized for <f> = — 7r/2 as 
^ 22/^31 < 0 [H|, and this implies that d e ff = 0 for 
the type II interaction oe —>■ e and the type 0 inter¬ 
action 00 —> o. For the cubic terms only included the 
XPM terms that involve the strong soliton and probe 
fields, and not the weak e-polarized SHG, SFG and DFG 
modes. We also note that even if BBO is anisotropic 
then the XPM term between two o-polarized modes is 

the same as the SPM term, i.e. Xxpm = Xspm = c n 
fl9| . In this identity Miller’s scaling is neglected, as it 
actually is in our code (a single frequency-independent 
Xj 3 J is used for each tensor element El)- However, for 
clarity let us keep them separated, as one could imagine 
cases where the probe is not the same polarization as 
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the soliton. For the probe we equivalently have 


id,E^ 




( W p) 


l/ o\ u 'p 

. p, 77 ,* iAfe° F ° £ 

+ ^- £? «.-a, p e “ s ~“ p 

i 7-1* 7-1 iAfc? p ° z 

+ 4.^ s +u p e 8 p 


+ 


3ca.e 


8 n 0 (u p )c 


Eu p [x'spmIAL p I + 2x XPM |£ l Ws | ]—0 


(A3) 


The phase-mismatch coefficients are 

Aklf G (9) = k e (2uj,9) - 2k 0 (uij) (A4) 

Afc w!+a; p ( 0 ) = kei^s +U p ,9) - k a (u s ) - k 0 (u p ) (A5) 
Afc°JiL(0) = k 0 {u s ) - k 0 (u p ) - k e (u! s - U) p , 0) (A6) 


where k e (u, 0) = n e (u, 9)u/c and as per usual n e (u, 9) = 
[cos 2 0/n 2 (w) + sin 2 0/n 2 (u;)] _1 / 2 where n e (w) is the 
BBO e-polarized refractive index. In what follows we 
drop the explicit dependence of the Afc’s and d e ff on 9 
and cj). 

Let us now write the basic plane-wave equations for 
the SHG, SFG and DFG e-polarized modes, disregarding 
irrelevant quadratic contributions (i.e. the ee —> e type 
0 interaction, which is very weak since it will be heavily 
phase mismatched and the relevant d e ff 0, Eq. (6)] is 
not as high as (I 33 in, e.g., LiNbOa), cubic nonlinearities 
(which we assume to a very good approximation to be 
irrelevant for the harmonics, as the phase mismatch is 
larger and thus their intensities too low) and chromatic 
dispersion (we are only interested in the nonlinear terms 
at the moment). 


id~E- 


2uijd e ff | 2 —'/A/;:' 


SHG 

Ef, e ‘“"'"J =0 


idzE^s —ujp ”h 


cn e (coj,0) 2 

{uj s - U p )d efi iAk°™„ p z = 0 

cn e {u s ~u p ,9) ^ 


(A7) 
= 0 
(AS) 


id z E, 


Z-^OJs+UJp 


(^S Wp)d e ff 77 77, 


cn e {u s + u p ,9) 


Eui.E u e 


+“» = 0 
(A9) 


With the usual cascading ansatz A kL 2n we can find 
the harmonic fields in the cascading limit [0 


E 2 u — ““ 


Ujd e ff 


cn e (uj,9)Ak^ G UJj 


„ 2 -iAk s * G z 

Ei.e 


(A10) 


Euj'.—uj^ — 


(Us Up') deft 


E t 




cn e (uj s - w p , 0 )AkSf? u 
(uj s ujp)d e ff 


7 -, 7 -,* iAk GFG z 

-E^E e --“p 


(All) 


cn e (uj s + 9)Ak^+ u 


77 . 771 — iAfc5 F 5„ z 

Eu.E^e 


(A12) 


When plugging these into Eqs. (lAB-dMl) we get 
3 u s „ r ( 3 ) 


id z E^j -f- 


8 n 0 {u s )c Euia [Aeff,SPM(^)l^ I ( Al 3) 

+ ^leff.XPMI Eu p |“] = 0 

id z E ulp + ■^-j^^E Up [x' e J SPM {u p )\E ulp \ 2 (A14) 

+ 2XgJ XPM |E Wa | ] = 0 


Here 


xSspm ( u j) =Xspm + X^ HG K) (A15) 

Xeff.XPM = XxPM A Xcasif G ( w s + w p) (A16) 
+ Xcasc DFG ( w s - w p) 


where 


(3),SHG/ \ _ 8 uJ j d ls 

Acasc ) ~ 


Xcasc SFG ( w s + W p ) = - 


3cn e (2uj, 9)Ak^ G 

3 (A17) 

_ 4(ta s +u p )d 2 eS 

3 cn e {u s + u p , 9)Ak^ G uip 


(A18) 


(3),DFG/ _ ) _ 4(cc., Up )d" t| 

x casc “ P )~ 3cne (u s - Up, 6)A kl :° F _% 

(A19) 

In the field normalized to the intensity case we have 
= E a , 3 .y / 2/e 0 cn o (w J ), j = s,p and 

id z A ulB + — A u1b [?r2^ I (w s )|A aJs |" (A20) 

+ 2< ff M KJ 2 ]=0 

*<4A, p + ^A Wp [nf™(w p )|A; p | 2 (A21) 

+ 2< ff M |A Wfl | 2 ]=0 

where 

«2™(Wj) =n 2™ rr ( a; j) + <£c(*i) (A22) 

< ff M =« r + + Wp ) + n» FG sc (u s - Wp ) 

(A23) 


and 


o (3) 

SPM / \ _ VXSPM 

Li.KerrWj 


4^o 


XPM 
^2,Kerr 


3x 


(3) 

XPM 


4ff 0 cn 0 (w s ) n D (wp) 


(A24) 

(A25) 
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Fig. 6. The induced cascaded nonlinear refractive indices 
vs. the SHG phase-mismatch parameter (controlled by the 
crystal angle 0), calculated for A s = 1.1 pm and \ v = 
1.65 pm. The domains are divided as follows: (l-|-4-|-5) Ef¬ 
fective self-focusing SPM regimes (nf™(w s ) > 0), so soli- 
tons cannot be excited; (2+3) Effective self-defocusing SPM 
regimes (nf™(o> s ) < 0), i.e. soliton regimes; (1+2+5) pos¬ 
itive XPM regimes (n 2 ™(w s ) > 0); (3+4) negative XPM 
regimes (nl+F (^s) < 0); (1+2) resonant nonlocal cascaded 
SHG regimes; (3+4+5) non-resonant nonlocal cascaded SHG 
regimes. 


while the cascading contributions are 

„shg / n = _ 2uJ o d l ff _ 

2 ’ cascl eo c?n2( Wj -)n e (2a; J -,6>)AfcSHG 

^2,easels T 8Jp) 

_ (++ + u P )d 2 eS _ 

E 0 c 2 n o (uj s )n o (u] p )n e (u] s + ui p , 9)Ak^ Up 

™2,cascK - Up) = 

_ (u s - u p )d 2 eff _ 

£ 0 c 2 n 0 (tu s )n 0 (uj p )n e (uj s - u p , 9)Ak^f^ p 

A requirement for the soliton to ” reflect” the probe is 
that the potential imposed by the soliton on the probe 
is scattering, i.e. a barrier. As the probe is launched 
in the anomalous dispersion regime, a scattering po¬ 
tential requires that the XPM effective nonlinear index 
of the probe (n*™) is negative, i.e. self-defocusing. As 
we showed above this has three contributions: = 

n +™r + «2?cS i c( w s+ w p) + n 2!cSc( w * viz - the ma¬ 

terial Kerr XPM as well as cascaded SFG and DFG be¬ 
tween the soliton and probe. 

In Fig. [5] we show how the crystal tuning angle 9 con¬ 
trols the various regimes; the self-defocusing soliton can 
only be excited in regimes (2+3). In turn the effective 
XPM term is only negative in regimes (3+4); at the 
boundary to regime (2) it becomes phase matched, and 
thus in (1+2) it is positive. Coincidentally the boundary 
between regimes (2) and (3) also marks the transition 
where the cascaded SHG becomes non-resonant, which 
it is in regime (3+4+5). In the non-resonant regime the 
cascading is ultrabroadband and induces minimal self- 
steepening on the soliton [3, IH, |25[ . 


(A26) 

(A27) 

(A28) 


Thus, for the wavelengths and tuning angles consid¬ 
ered here, ^2 easels — w p) is negligible due to a large 
phase mismatch. Therefore the XPM scattering poten¬ 
tial sign and magnitude is by and large a competition 
between the material Kerr XPM effect and the cascaded 
SFG effect. 
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